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PREFACE 


The  turbulence  data  presented  in  this  paper  were  collected  using 
sensors  known  as  "fat-wire  microtemperature  probes."  These  sensors  were 
a relatively  new  version  of  the  more  widely  used  "fine-wire"  probes. 

The  validity  of  fat-wire  probe  data  has  been  questioned  and  discussed 
since  a presentation  of  the  data  at  an  OSA  conference  on  turbulence  in 
August,  1977.  The  concern  centered  on  probability  distribution  functions 
calculated  from  the  data,  which  did  not  exhibit  the  expected  lognormal 
behavior.  Discussions  on  the  performance  of  fat-wire  probes  have  since 
taken  place  among  RADC,  NOAA  and  AFWL.  The  general  conclusion  was  that 
fat-wire  probe  data  are  valid  if  a correction  is  made  for  the  probe's 
limited  frequency  response.  In  this  report,  the  author  offers  three 
possible  explanations  for  the  discrepancy  between  the  theoretical  and 
experimental  probability  distributions. 

The  report  is  primarily  concerned  with  the  presentation  of  data  and 
certain  analytical  work.  Remarks  on  the  physical  meaning  of  the  data  are 
occasionally  made,  but  final  interpretation  is  left  for  the  reader  to 
make. 

The  assistance  provided  by  Major  John  Bradham  of  the  Staff  Meteoro- 
logy Office  at  RADC  is  appreciated.  Special  appreciation  is  expressed  to 
David  Youmans  for  his  efforts  which  kept  the  sensors  in  operation  through 
a cold  winter. 
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I.  INTRODUCTION 
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Atmospheric  degradation  of  an  optical  wavefront  is  an  important,  con- 
sideration to  researchers  concerned  with  the  optical  transmission  or 
reception  of  information.  The  surveillance  and  high  energy  laser  commu- 
nities are  now  developing  systems  aimed  at  minimizing  the  distortions 
produced  by  atmospheric  turbulence.  Turbulence  and  beam  degradation  data 
have  helped  formulate  and  refine  the  atmospheric  models  used  as  a basis 
for  the  competent  design  of  these  systems.  In  fact,  the  collection  of 
turbulence  data  during  optical  experiments  has  become  an  almost  vital 
requirement  at  test  ranges. 

Although  the  collection  of  turbulence  data  during  optical  experiments 
is  an  improvement  in  procedure  over  previous  years,  several  problems 
remain.  To  schedule  the  testing  of  an  instrument  under  specified  tur- 
bulence conditions,  the  researcher  has  had  to  guess  when  such  conditions 
would  occur.  His  guess  would  probably  be  based  on  personal  experience  and 
the  published  results  of  a few  short- term  experiments.  A need  has  there- 
fore existed  for  several  years  for  an  experiment  which  would  establish  a 
data  base  that  related  turbulence  strength  to  local  conditions,  such  as 
time  and  weather.  This  experiment  would  also  demonstrate  an  ability  to 
describe  "characteristic"  behavior  {or  behaviors)  of  atmospheric  turbulence 
near  the  ground,  in  spite  of  the  atmosphere's  inhomogeneous  nature  and 
statistical  nonstationarity. 


A long-term  experiment  measuring  atmospheric  turbulence  was  therefore 
conducted  at  the  Rome  Air  Development  Center  (RADC).  It  was  quickly 
established  that  the  "typical"  diurnal  turbulence  cycle  observed  by  others^ 
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actually  depended  on  the  local  weather,  In  particular  on  sky  conditions. 
The  experiment  had  to  be  terminated  after  four  months,  but  much  had 
already  been  learned.  Monthly  averages  of  the  diurnal  cycle,  probability 
distributions  and  other  statistics  have  been  evaluated  and  are  discussed 
in  this  report. 


II.  EXPERIMENTAL  APPROACH 


A.  LOCATION  AND  INSTRUMENTATION 

The  data  were  collected  at  the  Advanced  Optical  Test  Facility  (AOTF) 
of  RADC  in  Verona,  New  York.  Data  were  collected  for  the  months  of 
February  through  May,  1975.  The  physical  layout  of  the  site  is  shown  in 
Figure  1.  The  flat  terrain  is  covered  mostly  by  grassy  fields  which 
range  in  height  from  0.1  m in  the  winter  to  1.5  m in  the  summer.  Snow 
reached  its  greatest  depth  at  0.6  m in  February,  but  traces  of  snow  were 
recorded  into  April.  Cloudy  skies  dominated  the  winter  months,  which  was 
typical  for  the  region  (see  Table  1 for  climatic  probabilities  of  cloud 
cover  2).  Actual  weather  conditions  recorded  at  the  Griffiss  AFB  weather 
station  (approximately  25  km  east  of  the  AOTF)  during  the  months  of  the 
experiment  are  found  in  Tables  2-6  later  in  the  reporf. 

The  instrumentation  for  the  experiment  consisted  of  microtemperature 
probes,  wind  speed  and  wind  direction  sensors,  and  an  ambient  temperature 
system.  The  micro temperature  probes  were  fat-wire  versions  of  the  fine- 
wire  probe  described  in  detail  by  Greenwood  and  Youmans  The  tem- 

perature-sensing element  of  the  fat-wire  probe  was  a platinum  wire  12.7 
\x\\  in  diameter,  as  compared  to  the  2.5  ym  diameter  wire  of  the  fine-wire 
probe.  Although  the  system  response  of  the  two  versions  were  different, 
Greenwood  and  Spencer5  found  a relation  which  corrected  the  fat-wire  probe 
values  to  give  an  accurate  measurement  of  c£,  the  temperature  structure 
constant.  The  wind  speed  and  direction  systems  were  standard  products  of 
the  Climet  Instruments  Company.  The  wind  speed  sensor  had  a 0-5  volt 
analog  output,  corresponding  to  0-25  m/s  velocities.  The  wind  direction 
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sensor  was  a 0-5  volt,  360  degree  system.  The  ambient  temperature  system 
was  EG&G's  Model  110S-M  Temperature  and  Dew  Point  Measuring  Set. 

Three  microtemperature  probes  were  mounted  in  a triangular  configu- 
ration on  a small  tower  approximately  120  meters  SSE  of  the  AOTF  building 
(Fig.  2).  The  probes  were  mounted  in  a "standing"  position,  with  the 
wire  element  3 m above  ground  level  (see  Figure  3).  Wind  speed  and 
direction  sensors  were  placed  on  the  same  tower  at  a height  of  2.6  m. 

The  ambient  temperature  sensor  was  located  2 m above  the  ground  on 
another  tower  26  meters  SSE  of  the  AOTF  building. 

B.  DATA  ACQUISITION  AND  PROCESSING 

Data  from  the  six  instruments  on  the  towers  were  sequentially  sampled, 
digitized  and  recorded  by  a "datalogger"  at  a nominal  rate  of  one  cycle 
(six  data  channels)  per  second.  The  datalogger  was  a 12-bit  system  made 
by  Analog-Digital  Data  Systems,  Inc.  Six  to  eight  days  of  continuously 
sampled  data  were  stored  on  a standard  ten  inch  magnetic  tape  at  556  bits 
per  inch.  A new  file  was  marked  on  the  tape  daily  (except  weekends)  at 
1200  EST.  All  magnetic  tapes  were  processed  at  RADC  on  a Honeywell  635/ 

645  computer.  A modified  version  of  a FORTRAN  program  written  by  Pattern 
Analysis  and  Recognition  Corporation6  was  used  to  reduce  the  raw  data  to 
30-minute  averages.  These  averages  included  the  mean  and  standard  devi- 
ation for  wind  speed,  wind  direction  and  ambient  temperature,  as  well  as 
c|  for  each  of  the  three  possible  probe  pairings.  The  averages  formed  the 
basis  for  almost  all  subsequent  analysis.  Appendix  A examines  the 
"smoothing"  effect  of  using  30-minute  periods  as  compared  to  shorter  5- 
minute  periods. 
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III.  EXPERIMENTAL  FINDINGS 


A.  EARLY  RESULTS  AND  PROBLEMS 

The  temperature  structure  constant  is  calculated  from  two  microtem- 
perature probes  from  the  relation^ 

2 


Cl  * fri2  + - 2ClJ]  rij‘  ' , V<r1j«Lo  (1) 

p p 

where  <jj  and  aj  are  the  variances  of  the  probes  with  separation  r^j, 

Cjj  is  the  probes'  covariance,  £0  is  the  "inner  scale"  of  Kolmogorov's 
inertial  subrange  (of  the  temperature  temporal  power  spectrum)  and  LQ  is 
the  "outer  scale."  The  inner  scale  is  typically  a few  millimeters  in 
size,  while  the  outer  scale  is  on  the  order  of  height  above  ground.  An 
analysis  shows  that  even  if  the  two  probes  are  not  accurately  calibrated, 
the  calculated  value  of  is  not  affected  significantly  (see  Appendix  B), 
Equation  (1)  applies  to  fine-wire  probes  - it  does  not  include  the 
frequency  response  correction  for  fat-wire  probes.  Multiplying  the 
results  of  equation  (1)  by 

1.2  (v/r-j j)° * 182  , v/rj j<20  (2) 

where  v is  the  wind  speed,  yields  the  "corrected"  value  for  from  fat- 
wire  probes.  The  correction  factor  given  by  (2)  underestimates  the 
needed  correction  for  values  of  v/r> 20.  At  v/r=20,  the  error  in  (2)  is 
ZJ%.  The  error  increases  to  9%  at  v/r=40.  Since  the  values  of  v/r 
which  occurred  during  the  experiment  were  rarely  greater  than  40,  the 
correction  given  by  (2)  is  considered  adequate. 

Analysis  of  initial  data  showed  that  the  smoothest,  most  regular 
diurnal  profiles  of  were  obtained  by  selecting  the  values  from  the 
probe  pairings  whose  line  of  separation  was  most  perpendicular  to  the 
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wind  direction.  This  apparent  dependence  on  wind/probe  orientation  has 
been  a topic  of  discussion  for  many  years,  but  no  experimental  studies 
have  been  reported.  The  probes  were  mounted  in  a triangular  configuration 
in  this  experiment  specifically  to  avoid  any  problems  with  wind/probe-pair 
orientation.  For  each  30-minute  average,  the  computer  program  which  re- 
duced the  raw  data  automatically  picked  the  "best"  pair  of  probes  according 
to  wind  direction. 

The  triangular  configuration  did  not  resolve  other  problems,  however. 
Spurious  effects  were  observed  when  the  wind  was  not  "good"  for  any  of 
the  three  probe  pairings  - specifically,  problems  occurred  when  the  wind 
direction  was  within  15°  of  the  angle  bisector  of  either  acute  angle  of 
our  (obtuse)  triangle  configuration.  Similar  effects  occurred  when  the 
winds  were  light  and  variable.  This  condition  was  identified  by  a 
standard  deviation  for  wind  direction  greater  than  30  degrees.  To  avoid 
the  apparent  wake  effect  problems  of  these  two  conditions,  was  cal- 
culated from  the  single,  windward  probe  using  the  relation'7 

c$  * 1.910V V 2/3  (3) 

Approximately  1 0-30%  of  the  48  averages  of  each  day  were  calculated 
using  equation  (3).  Since  the  outer  scale,  Le,  is  not  a static  quantity, 
the  use  of  equation  (3)  undoubtedly  introduced  an  error  in  the  calculated 
value  of  C(j.  Ideally,  t0  would  have  been  derived  from  data  immediately 
preceding  the  period  when  equation  (3)  was  needed.  This  proved  to  be 
impossible  for  such  a long-term  experiment.  Average  values  of  were 
calculated  instead  by  substituting  60  known  values  of  c£  from  various 
times  of  day  into  equation  (3).  The  mean  values  obtained  in  this  manner 
were 
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Lq  = 5.42  m (night) 

L0  = 1 .26  m (day) 

which  only  apply  tc  fat-wire  probes.  Comparison  of  known  values  with 
values  calculated  from  equation  (3)  showed  that  the  error  was  usually 
less  than  a factor  of  2.  Also,  each  daily  turbulence  profile  was  examined 
to  insure  that  the  single-probe  calculations  remained  reasonable. 

Another  early  observation  was  the  dependence  of  the  daily  turbulence 
profile  on  sky  conditions.  The  turbulence  profile  for  each  24  hours  of 
data  was  eventually  classified  as  belonging  to  three  categories.  Type  I 
turbulence,  as  defined  here,  exhibits  a clear  diurnal  cycle.  Figure  4, 
a typical  example,  shows  the  turbulence  levels  for  the  24-hour  period 
beginning  at  noon  (local  standard  time)  on  1 February  1975  (note  that 
all  similar  references  to  specific  dates  also  imply  noon  as  the  start  of 
the  period). 

For  Type  I turbulence,  the  profile  peaks  at  midday,  but  it  rapidly 
falls  off  just  prior  to  sunset.  After  this  dip,  it  usually  regains  much 
of  its  strength,  but  subsequently  decreases  during  the  night.  A few 
nights  show  a moderate  peak  later  similar  to  the  daytime  peak. 

Type  I turbulence  is  usually  associated  with  clear  skies,  scattered 
clouds  or  perhaps  a light  overcast.  Type  II  turbulence  is  a degraded 
form  of  Type  I turbulence.  Type  II  turbulence  profiles  are  not  as  smooth, 
with  more  spikes  and  dips.  The  profile  for  29  April  1975  (Figure  5)  is  an 
example.  Strong  midday  turbulence  still  occurs,  and  the  "sunset"  cip 
can  still  be  identified.  This  type  of  turbulence  is  generally  associated 
with  sky  conditions  ranging  from  partly  cloudy  to  moderate  overcast. 


Type  III  turbulence  profiles  are  even  further  degraded.  As  seen  in 


Figure  6,  the  characteristic  profile  is  difficult  to  identify.  This  type 


of  turbulence  is  usually  associated  with  total  cloud  cover  and  possibly 


some  precipitation.  The  occurrence  of  precipitation  makes  the  interpre- 


tation of  Type  III  turbulence  a difficult  proposition,  if  the  values  have 


any  meaning  at  all . 


Examination  of  Figures  4,  5 and  6 reveal  the  characteristic  features 


of  the  different  turbulence  categories.  Strong  daytime  turbulence 


dominates  the  profiles  for  both  1 Feb  75  and  29  Apr  75.  Pronounced 


"sunset"  and  "sunrise"  dips  are  also  evident  in  both.  The  sunset  dips 


are  30  to  60  minutes  before  actual  sunset  (note  that  the  data  points  are 


plotted  against  the  beginning  time  of  each  30-minute  averaging  period). 


The  sunrise  dips  occur  near  or  after  the  actual  time  of  sunrise.  Whereas 


the  sunset  dip  is  universal  to  all  Type  I and  Type  II  turbulence  profiles, 


the  sunrise  dip  is  not  as  common.  This  will  be  evident  in  the  monthly 


averages  discussed  later.  The  profile  of  29  Apr  75  has  been  classified 


as  Type  II  because  of  its  jumpy  behavior  between  1200  and  2000  EST,  not 


because  of  the  strong  nighttime  turbulence.  Such  nighttime  turbulence 


was  part  of  a trend  which  will  be  discussed  with  the  monthly  averages. 


Note  that  about  0600,  the  profile  is  observed  transitioning  to  Type  I 


again.  The  profile  from  20  Feb  75  demonstrates  the  irregular  behavior 


of  Type  III  turbulence,  although  sunset  and  sunrise  dips  might  still  be 


identified. 


Certain  aspects  of  the  experiment  should  be  clarified  before  proceding 


further.  The  classification  of  the  turbulence  "type"  resulted  from  a 


series  of  subjective  decisions  made  over  an  extended  period  time. 
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No  single  criterion  was  used  to  distinguish  between  types.  The  decision 
to  use  three  turbulence  categories  was  based  on  the  desire  to  obtain 
better  statistics  than  those  which  could  be  obtained  from  a simpler  good/ 
bad  classification.  Periods  of  low  wind  speeds  (below  1 m/s)  were  not 
omitted,  in  contrast  with  short-term  experiments  by  others.  Light 
precipitation  sometimes  occurred  even  in  connection  with  Type  I turbulence 
profiles.  The  actual  weather  conditions  recorded  at  Griffiss  AFB  Weather 
Station  are  summarized  in  Tables  2-6,  according  to  month,  day,  relative 
time,  and  turbulence  classification. 

B.  MONTHLY  RESULTS 

The  monthly  results  were  computed  two  different  ways.  In  both  cases, 
the  turbulence  values  from  the  30-minute  blocks  of  raw  data  were  combined 
into  60-minute  blocks  corresponding  to  the  hours  of  the  day.  The  initial 
approach  was  to  calculate  the  arithmetic  mean  ard  the  standard  deviation 
for  each  period.  The  second  approach  was  to  obtain  a median  value  for 
each  60-minute  period.  The  median  was  obtained  from  a discrete  cumulative 
probability  distribution  which  was  calculated  for  each  period.  The  number 
used  here  as  the  "median"  is  the  discrete  value  below  which  lies  50  percent 
of  the  distribution. 

The  monthly  averages  of  the  daily  turbulence  profiles  were  calculated 
and  plotted  in  terms  of  the  temperature  structure  constant,  Cfj.  The 
refractive  index  structure  constant,  C^,  can  be  found  from  these  values 
using  a simple  relation  (see  Greenwood  and  Youmans,  p.  55).  Wind  speed 
and  ambient  temperature  associated  with  periods  of  Type  I turbulence  were 
also  averaged. 

Several  important  findings  are  evident  in  the  monthly  averages  of  Type 
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I turbulence  {see  Figures  7-14).  The  strength  of  the  turbulence  during 
the  daylight  hours  Is,  of  course,  the  outstanding  characteristic  of  the 
profiles.  As  the  hours  of  sunshine  increased  from  February  to  May,  the 
strong  daytime  turbulence  also  lengthened.  The  strongest  values  always 
occurred  between  1100  and  1300  EST.  The  peak  values  of  the  four  monthly 
profiles  are  almost  constant  even  though  ambient  temperatures  change 
dramatically.  The  midday  peaks  of  c£  vary  only  between  0.48  °C2m"2/3 
for  March  and  0.37  °C2nf2/3  for  February.  Figures  11-14  show  that  the 
average  temperatures  were  over  20°C  warmer  in  May  than  they  were  in  Feb- 
ruary. 

Other  characteristics  and  trends  can  be  found  in  the  plots  of  Type  I 
turbulence.  Nighttime  turbulence  drops  to  its  lowest  strengths  for 
February  and  March,  but  it  is  stronger  for  the  warmer  months  (see  Figure 
15).  February  has  the  lowest  nighttime  value  of  c£  with  0.004  °C2m-2/3, 
while  May's  lowest  value  at  night  (excluding  the  "sunset  dip")  is  only 
down  to  0.023  °C2m"2/3.  If  this  trend  continues  into  the  summer  months, 
it  would  support  the  common  belief  that  the  best  nights  for  "seeing"  are 
cold  winter  nights.  The  sunset  dip  seen  in  the  profiles  occurs  at  later 
times  as  local  sunset  gets  later.  The  value  of  this  dip  is  fairly  constant 
from  month  to  month,  varying  only  between  0.0085  and  0.017  0C2m"2^3. 

There  is  no  "sunrise  dip"  which  can  be  clearly  identified  in  all  the 
monthly  profiles,  although  It  does  occur  occasionally  on  certain  days 
(see  Figures  4 and  5).  The  strongest  sunrise  dip  in  a monthly  average 
occurred  in  April  for  the  0600-0700  time  period  (see  Figure  9),  which  is 
roughly  30  minutes  after  the  median  time  of  sunrise  for  that  data. 
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Upon  examining  the  averages  of  Type  II  turbulence  (Figures  16-19), 
slightly  "rougher"  profiles  are  seen.  February's  profile  is  omitted  from 
the  following  discussion  and  is  addressed  in  the  next  paragraph.  As 
compared  to  the  Type  I profile,  the  Type  II  profile  for  March  is  different 
by  a later  build-up  and  an  earlier  fall  off  in  the  daytime  turbulence.  The 
April  and  May  Type  II  profiles  differ  from  Type  I by  irregular  jumps  in 
the  turbulence  values  during  daylight  hours.  All  three  months  show  stronger 
turbulence  at  night  for  Type  II  profiles  than  for  Type  I profiles,  while 
maintaining  the  trend  toward  stronger  nighttime  turbulence  in  the  warmer 
months  (see  Figure  20). 

The  February  profile  which  has  been  plotted  as  Type  II  turbulence 
(Figure  16)  obviously  does  not  fit  the  usual  definition.  Only  three 
unusual  days  in  February  were  classified  as  Type  II,  so  statistics  are 
extremely  poor.  Figure  16  has  been  Included  only  for  completeness.  It 
is  strongly  recommended  that  no  conclusions  should  be  drawn  from  any 
reference  to  the  February  Type  II  turbulence  data. 

When  all  the  turbulence  types  are  combined  for  an  "all  weather"  average, 
an  obvious  effect  takes  place.  The  diurnal  cycle  of  the  Type  I turbulence 
profile  is  still  exhibited,  but  the  midday  peak  is  weaker  and  the  nighttime 
turbulence  is  stronger  (see  Figures  21  and  24).  This  suggests  that  tur- 
bulence during  cloudy  weather  changes  less  from  night  to  day.  April's 
"all  weather"  profile  shows  the  least  change  from  the  Type  I profile  be- 
cause there  were  no  Type  III  turbulence  data  available  for  that  month. 

That  lack  was  probably  due  to  the  heavy  rain  showers  in  April,  which  would 
have  broken  the  micro temperature  probes'  sensing  elements  more  readily 
than  would  February's  and  March's  snows.  The  drastic  degradation  in 
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February's  profile  Is  due  mostly  to  Type  III  turbulence  (snow?)  and  only 
slightly  to  the  Type  II  turbulence  profile. 

The  standard  deviations  of  the  monthly  averages  of  turbulence  are 
plotted  as  short  dashes  in  many  of  the  figures.  The  standard  deviations 
are  often  as  high  as  the  average,  which  is  to  be  expected  from  the  non- 
stationarity  of  the  atmosphere.  For  Type  I turbulence  profiles,  the 
standard  deviations  tend  to  be  the  smallest  (relative  to  the  mean)  during 
the  day.  Almost  all  their  values  at  night  are  on  the  same  order  as  the 
means.  The  situation  almost  reverses  for  Type  II  turbulence,  with  the 
smallest  standard  deviations  occurring  late  at  night.  All  the  standard 
deviations  of  the  "all  weather"  profiles  are  high.  In  the  plots,  the 
smaller  standard  deviations  indicate  that  similar  values  of  turbulence 
were  observed  from  day  to  day,  in  spite  of  the  nonstationarity  of  the 
atmosphere.  The  large  standard  deviations  introduce  ar.  uncertainty  to 
the  average  values  which  must  not  be  ignored.  This  does  not  alter  the 
general  interpretation  of  the  monthly  averages,  but  rather  it  reinforces 
the  appreciation  for  the  statistical  nature  of  this  type  of  data. 

Because  of  the  large  standard  deviations  associated  with  the  monthly 
averages,  the  median  values  were  also  examined.  The  monthly  median  values 
of  Type  I turbulence  are  plotted  in  Figures  25  - 28.  They  show  that 
the  medians  for  4-6  hours  around  midday  were  almost  the  same  as  the 
arithmetic  means,  but  the  nighttime  median  values  were  roughly  2-5  times 
lower  than  the  corresponding  mean  values.  This  was  expected  from  the 
behavior  of  the  standard  deviations.  For  the  other  types  of  turbulence, 
the  use  of  median  values  instead  of  average  values  caused  more  noticeable 
changes  in  the  shapes  of  the  profiles. 
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IV.  STATISTICAL  ANALYSIS 


A.  PROBABILITY 

The  statistical  distribution  of  c£  values  has  been  debated  for  several 
years.  Most  of  these  discussions  have  been  based  on  a few  days  of  data, 
whereas  four  months  of  data  were  collected  for  this  experiment.  It  was 
hoped  that  the  question  of  distributions  could  finally  be  answered. 

An  initial  look  at  the  data  was  made  for  each  of  the  three  types  of  tur- 
bulence for  day  versus  night.  The  distributions  were  not  well-defined, 
so  it  was  suggested®  that  data  collected  around  dawn  and  dusk  be  excluded. 

The  final  analysis  calculated  the  discrete  probability  distribution 
function  and  the  cumulative  probability  function  from  all  the  30-minute 
c£  values  according  to  turbulence  type  and  time  period.  The  day  was 
divided  into  the  following  four  time  periods: 

night  - 1900  to  0400  EST 

day  - 0800  to  1530  EST 

dusk  - 1530  to  1900  EST 

dawn  - 0400  to  0800  EST 

The  long  times  allowed  for  "dusk"  and  "dawn"  were  specified  to  prevent 
contamination  of  the  day  and  night  periods  by  the  sunset  and  sunrise 
dips.  The  plots  of  the  probability  functions  for  Type  I turbulence.  Type 
II  turbulence  and  "all  weather"  conditions  are  included  in  Appendix  C. 

After  examining  some  of  the  calculated  functions,  it  is  evident  that 
there  is  no  distinct  lognormal  distribution  (since  log  is  plotted,  a 
Gaussian  distribution  is  the  actual  focus  of  the  search).  The  daytime 
category  has  a probability  distribution  which  appears  to  be  characteristic 
for  most  months.  The  distribution  is  usually  strongly  skewed  to  the  left, 
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and  it  sometimes  has  a Gaussian-1 ike  peak.  An  extreme  example  of  this  is 
the  daytima  distribution  for  Type  I turbulence  in  May  (see  Figure  29). 

The  distribution  for  the  nighttime  category  are  less  consistent,  except 
that  they  have  a much  milder  skew  than  the  daytime  distribution.  For  all 
types  of  turbulence,  most  dusk  and  dawn  categories  have  wide,  flat  distri- 
butions. 

The  quantitative  statistics  associated  with  the  distributions  are 
summarized  in  Tables  7-9.  A true  Gaussian  distribution  would  have  a 
skewness  of  0 and  a kurtosis  of  3.  Very  few  of  the  distributions  are 
quantitatively  close  to  a Gaussian  form.  The  daytime  distributions  for 
Type  I turbulence  are  particularly  aberrant.  This  was  not  expected.  From 
the  tables,  the  distributions  of  log  C?  values  which  come  closest  to 
Gaussian  statistics  are  February's  daytime  Type  II  turbulence  (described 
earlier  as  atypical  and  not  usable)  and  March's  nighttime  "all  weather" 
turbulence.  The  March  data  is  plotted  on  a probability  scale  in  Figure 
30.  A Gaussian  distribution  plotted  on  this  scale  would  give  a straight 
line.  For  a different  perspective,  the  March  data  is  plotted  conventionally 
with  a Gaussian  fit  (Figure  31).  This  example  is  not  bad  by  itself,  but 
unfortunately  it  is  the  curve  closest  to  Gaussian  form  from  all  types  of 
turbulence. 

Skimming  through  all  the  distributions  plotted  in  Appendix  C and  not 
consulting  the  statistics  tables,  the  reader  might  want  to  pick  the 
distribution  for  May  Type  II  turbulence  at  night  as  closest  to  a Gaussian 
distribution.  A plot  of  it  in  Figure  32  confirms  that  this  distribution 
was  worse  than  the  one  picked  from  the  tables. 

To  examine  the  difference  between  the  observed  distributions  and  a 
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true  Gaussian  distribution,  mean  squared  errors  were  calculated  for  the 
coefficients  of  skewness  and  kurtosis.  Table  10  summarizes  the  results. 

The  figures  in  Table  10  suggest  that,  for  overall  values,  the  "all  weather" 
distributions  are  the  closest  to  Gaussian  form,  while  the  Type  I turbulence 
distributions  are  the  least  Gaussian.  After  breaking  the  results  down 
according  to  time,  however,  the  apparent  differences  fade.  The  dusk  and 
dawn  values  will  continue  to  be  ignored,  since  their  characteristic  dips 
contaminate  the  distributions.  Type  I turbulence  for  daylight  hours  is 
definitely  the  least  Gaussian-like  of  the  daytime  data.  The  mean  squared 
difference  in  skewness  is  1.50,  which  is  due  to  consistently  strong  skews 
to  the  left.  The  high  error  in  kurtosis  is  almost  entirely  due  to  strongly 
peaked  distributions  for  February  and  May.  For  the  daytime  periods,  Type 
II  turbulence  distributions  have  the  lowest  skewness,  while  "all  weather" 
conditions  have  the  lowest  error  in  kurtosis.  At  night,  Type  I turbulence 
distributions  are  the  least  skewed  but  the  worst  in  kurtosis.  The  least 
error  in  kurtosis  at  night  is  found  for  the  "all  weather"  conditions. 

By  now  it  is  clearly  obvious  that  the  predicted  laynormal  behavior  is 
not  confirmed  by  this  experiment.  The  most  Gaussian-like  distributions 
were  expected  from  the  daytime  values  for  Type  I turbulence,  but  these 
distributions  are  strongly  skewed.  Since  they  show  an  unexpected,  abrupt 
cutoff  of  high  values,  the  problem  appears  to  be  related  to  the  data  rather 
than  theory.  There  are  at  least  three  possible  explanations  for  such 
behavior.  The  first  possibility  is  that  30-minute  averaging  periods  are 
too  long  and  the  data  is  smoothed  too  much  to  observe  the  true  distribution. 
The  implication  is  that  periods  of  high  turbulence  are  short-lived  or  at 
least  intermittent.  Another  explanation  is  that  sampling  the  raw  signal 
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at  a rate  of  1 Hz  is  too  slow  to  measure  high  turbulence.  A third 
possibility  contends  that  the  limited  frequency  response  of  fat-wire 
probes  (approximately  0-20  Hz)  simply  prevents  accurate  measurements  of 
strong  turbulence. 

Only  detailed  experimental  data  can  resolve  the  issue.  Results  from 
Appendix  A confirm  that  averaging  time  does  affect  the  shape  of  the 
distributions.  By  shortening  the  averaging  times  to  5 minutes,  the 
abrupt  cutoff  of  the  distribution  lengthens  out  into  a "tail",  thereby 
improving  distribution  shape.  Sampling  rate  would  become  more  important 
as  averaging  periods  were  shortened  to  one  minute  or  less.  In  1977  Norton 
and  Walters^  reported  reasonable  agreement  between  predicted  and  measured 
probability  distributions.  Their  data,  also  collected  with  fat-wire  probes, 
calculated  turbulence  strength  for  10-second  averages  over  15  minute 
windows  in  time.  These  positive  results  further  support  averaging  time 
and  sampling  rate  as  major  factors  in  the  analysis.  The  data  needed  to 
evaluate  the  performance  of  fat-wire  probes  in  strong  turbulence  is  now 
being  collected  at  the  AOTF.  Atmospheric  nonstationarity,  which  is  another 
consideration,  is  apparently  not  the  problem  because  it  would  smear  out  the 
probability  distributions.  That  would  be  exactly  the  opposite  of  the 
effect  observed.  The  eventual  answer  might  be  that  the  three  factors  - 
averaging  time,  sampling  rate  and  frequency  response  of  fat-wire  probes  - 
can  each  limit  the  measurement  of  strong  turbulence  and  affect  the  shape 
of  probability  distributions. 
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B.  POWER  SPECTRA 


The  turbulence  values  were  next  analyzed  for  any  possible  hidden  cyclic 
behavior,  as  well  as  for  different  power  spectral  forms  according  to  type 
of  turbulence.  The  values  of  the  power  spectrum  and  autocovariance  function 
were  calculated  for  35  different  daily  profiles,  each  consisting  of  a full 
24  hours  of  data.  These  were  divided  according  to  month  and  type  of 
turbulence  to  identify  any  common  features. 

The  power  spectrum  and  the  autocovariance  values  were  calculated  in 
' the  following  way  11.  (i)  a 48-point  discrete  Fourier  transform  was 

i 

j applied  to  the  input  data  (with  the  average  already  removed);  (2)  a 

Hamming  window  was  applied  as  a 3-point  convolver;  (3)  the  raw  power 
| spectrum  was  taken  to  be  the  "modulus  squared"  of  the  resulting  sequence; 

l 

i (4)  the  raw  spectrum  was  averaged  over  adjacent  bands  (values)  two  wide; 

(5)  the  spectrum  was  normalized  such  that  the  total  power  equalled  the 
variance;  and  (6)  the  autocovariance  was  found  as  the  inverse  Fourier 
transform  of  the  power  spectrum.  Some  of  these  results  are  plotted  in 
Fiqu>es  33-45. 

Common  features  are  evident  ir  the  plots.  Examining  the  power  spectra 
first,  we  see  that  Type  I turbulence  spectra  have  the  sharpest  and  strongest 
central  peaks.  There  is  usually  at  least  one  small  peak  in  the  middle 
frequencies  at  about  0.5  hr“^.  For  Type  II  turbulence  the  central  peak  is 
wider  with  somewhat  less  power,  whereas  the  higher  frequencies  have  more 
power.  For  the  spectra  of  Type  III  turbulence,  the  central  peak  is  also 
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wide,  with  more  energy  shifted  out  to  the  higher  frequencies.  These 
features  of  the  power  spectra  indicate  that,  for  good  weather  with 
corresponding  "good"  turbulence  profiles,  the  values  are  dominated  by 
a low  frequency  component.  This  component  has  a period  of  approximately 
12  hours,  which  (conveniently)  corresponds  to  a day/night  pattern.  With 
cloudier  weather,  the  low  frequency  component  becomes  less  dominant  while 
the  higher  frequencies  increase  in  power.  The  reduction  in  low  frequencies 
might  reflect  a lessened  change  between  day  and  night  turbulence.  Relying 
solely  on  power  spectra  of  sequences  of  48  30-minute  c|  values,  though, 
it  is  difficult  to  differentiate  the  types  of  turbulence. 

Differences  are  more  evident  in  the  autocovariance  function.  Examining 
the  shape  of  the  autocovariance  functions,  we  see  that  Type  I turbulence  is 
clearly  different  from  the  others.  For  Type  I turbulence,  the  function  has 
a form  intermediate  between  a Gaussian  and  a triangle  function.  Relative 
to  other  types  of  turbulence  for  the  same  month,  Type  I turbulence  also 
tends  to  have  stronger  autocovariance  values.  Type  II  and  Type  III 
turbulence  exhibit  a flatter  form  for  the  function,  although  some  of  the 
Type  III  plots  display  a distinctive  shape  resembling  the  sine2  function. 
The  conclusion  is  that  the  regular,  smooth  nature  of  Type  I turbulence 
profiles  is  a result  of  a definite,  positive  correlation  between 
successive  values12.  Such  correlation  either  deteriorated  or  never 
existed  for  the  other  turbulence  types. 


V.  SUMMARY 


The  data  has  now  been  presented  and  briefly  discussed.  The  results, 
which  are  specifically  applicable  to  work  at  the  AOTF,  should  be  useful 
at  other  locations.  Yet,  what  conclusions  can  be  drawn  from  all  this 
data? 

First,  the  scheme  for  classifying  turbulence  into  three  basic  categories 
appears  reasonable.  The  three  types  of  turbulence  are  associated  with 
weather  conditions  which  could  almost  be  termed  "good",  "fair"  or  "poor". 

The  standard  deviations  of  the  turbulence  values  for  monthly  results  were 
highest  for  the  "all  weather"  combination,  still  high  for  the  average 
daytime  values  of  Type  II  turbulence,  but  relatively  low  for  daytime  values 
of  Type  I turbulence.  The  need  for  an  intermediate  turbulence  form  (Type 
II)  is  supported  by  the  analysis  of  autocovariance  functions,  which  showed 
a discernible  difference  between  Type  I and  Type  II  turbulence  in  spite 
of  similarities  in  their  profiles. 

Type  I turbulence  exhibits  some  exceptional  features  in  the  monthly 
profiles.  The  midday  peak  is  roughly  the  same  over  four  months.  The 
sunset  dip  also  remains  fairly  constant.  In  contrast,  turbulence  strength 
at  night  increases  with  the  warmer  months.  These  findings  should  be 
interesting  to  investigators  studying  turbulence  from  a heat  budget 
standpoint. 

Although  the  probability  distributions  do  net  display  the  predicted 
lognormal  shape,  reasons  are  offered  which  might  explain  the  discrepancy. 

The  experimental  curves  are  valid  within  the  limitations  apparently  imposed 
by  the  manner  in  which  the  data  were  collected  and  reduced. 

Other  conclusions  might  be  drawn  from  the  data,  but  they  are  left  to 


the  reader.  Hopefully,  this  report  will  not  mark  the  only  use  of  the  data. 
An  attempt  has  been  made  to  present  enough  data  to  allow  immediate  use  by 
others  who  could  not  otherwise  Invest  the  required  time  and  resources. 


t==5=^S=)  f===) 


- 7m  TREES 


5 ■ 

n*  V 


BROOK 


HAY  FIELDS 


jjj  INSTRUMENT 
~T“  .STATIONS' 


WP? 

fe  D 

2 " 


O.l-  Im 
GRASS 


DIRT  ROAD 

I \ 


~9m 

TREES 


•*.,  ~5m 
:.  V TREES, 

sparse 


u/ . 


~7in 

TREES 


• I m 
BRUSH 


AOTF— i ^ BRUSH 

tO  '•  • 


-1  m 
TREES 


— 2 m 
BRUSH 


•10m  TREE  ROW 


7m  TREES, 
SPARSE 


~l/2  m 
LOW  GRASS 


Figure  1.  Terrain  around  the  AOTF.  Instrument  towers;# 
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Figure  2.  Instrument  tower 


Example  of  Type  III  turbulence  (20  Eeb  75). 


Figure  9.  Mean  (+)  and  std.  dev. 
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Figure  12.  Average  wind  speed  and  ambient  temperature  recorded 
during  periods  of  Type  I turbulence:  March  1975. 
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Figure  13 
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. Average  wind  speed  and  ambient  temperature  recorded 
during  periods  of  Type  I turbulence:  April  1973. 
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Fig-are  15.  Mean  values  of  -Type  I turbulence:  February  - May  1975 


Figure  16.  Kean  (+)  and  std.  dev.  (-)  for  Type  II  turbulence:  February  1975 


?igure  17.  Kean  (+)  and  std.  dev.  (-)  for  Type  II  turbulence:  March.  1975 


Figure  19.  Mean  (+)  and  std.  dev.  (-)  for  Type  II  turbulence:  May  1975 


Figure  20.  Mean  values  of  Type  II  ‘turbulence:  March.  - May  1975 


Figure  21.  Mean  (+)  and  std.  dev.  (-)  for  "all  weather" : February  1975 


May  1975 


Figure  25.  Median  values  for  Type  I turbulence:  February 


Figure  26.  Median  values  for  Type  I turbulence:  March. 


Figure  30.  Cumulative  probability  function  for  "all  weather":  March 
night  values. 


Figure  31.  “All  weather”  turbulence  values  (March  1975,  night  values)  compared  to  its 
Gaussian  fit  (solid  line). 
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Figure  55.  Power  spectrum  (+)  and  auto covariance  function  (-) 
for  11  Feb  75  (Type  III,  typical). 
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Figure  37.  Power  spectrum  (+)  and  auto covariance  function  (-) 
for  8 Mar  75  (Type  I). 
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Figure  38.  Power  spectrum  (+)  and  autocovariance  function  (-) 
for  6 Mar  75  (Type  II). 
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Figure  40.  Power  spectrum  (+)  and  autocovariance  function  (-) 
for  10  Apr  75  (Type  I). 
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Figure  42.  Power  spectrum  (+)  and  auto covariance  function  (-) 
for  2 May  75  ( Type  I). 
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Figure  43.  Power  spectrum  (+)  and  auto oo variance  function  (-) 
for  4 May  75  (Type  I). 
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TABLE  1.  CFLOS  PROBABILITY  {%) 


X i 


1 


i 


* i 


Time  (LST) 


Month 

03-05 

09-11 

15-17 

21-23 

Jan 

33 

35 

35 

35 

Feb 

38 

39 

38 

39 

Mar 

42 

40 

58 

44 

Apr 

46 

40 

41 

48 

May 

31 

45 

46 

54 

Jun 

56 

51 

51 

61 

Jul 

61 

54 

55 

65 

Au  g 

59 

52 

56 

65 

Sep 

60 

52 

55 

62 

Oct 

59 

50 

53 

59 

Nov 

38 

34 

55 

38 

Dec 

33 

32 

34 

34 

t 


Note:  Viewing  angle  is  90°  elevation. 
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TABLE  2.  LEGEND  EOR  WEATHER  CONDITIONS 


0 Clear  skies 
8 Scattered  clouds  ( <50$) 

0 Broken  cloud  cover  (50-99$) 

@ Overcast  (100$) 

WX  Cloud  cover  not  observable 
R Rain 
S Snow 
L Drizzle 
E Eog 
H Haze 
K Smoke 
T Thunder 
Z_  Freezing  ... 

_W  ...  Shower 
B_  Blowing  . . . 

+ (heavy) 

(light) 
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TABLE  3.  WEATHER  - FEBRUARY  1975 


Late  Sky 

(Type  I Turbulence) 

31  Jan  0 to  9,  @ 

1 Feb  0,  0 

2 8 

3 0;  later  0 

(Type  II  Turbulence) 

6 Feb  © 

7 @ 

8 @ 

18  @ 

("Quiet  Turbulence") 

4-  Feb  © 

5 © 

(Type  III  Turbulence) 


Precipitation 
some  SW- 
some  SW- 


frequent  S-,  F 
frequent  SW-,  BS 
SW- , S-,  BS 
F,  R-,  L-,  SW- 

late  SW-,  S- 
S-f  F 


Snow  Depth 


<1  cm 
<1 
<1 
<1 

22-25  cm 
25 
25-27 
15-10 

<1-5  cm 
12-21 


8 Feb 

© 

late  SW- 

25-27 

9 

@ to  0 to  0 

S-,  SW- 

27-30 

10 

@ 

SW-,  s- 

30-33 

11 

0,  © to  0 to  © 

S-,  SW-,  F 

33-35 

12 

e,  wx 

late  SW- , SW+ , BS , SW 

35-50 

13 

wx,  ©,  © 

SW- 

50-55 

19 

© 

SW- 

10 

20 

@ 

10 

26 

mostly  ® 

SW- 

<1 

27 

mostly  © 

SW- 

<1 

28 

mostly  @ 

SW-,  T,  SW 

<1 

TABLE  5.  WEATHER  - APRIL  1975 


Date 

Sky 

(Type  I Turbulence) 

8 Apr 

@ to  0 to  0 

S 

8,0  to  © , @ 

10 

@ to  0,  0 

11 

0,  0 

14 

8 to  @ 

15 

@ 

16 

©to  0 

21 

0 to  0 to  0 

22 

0 to  0 

25 

0 to  @ 

28 

0,  @ 

(Type  II  Turbulence) 

1 Apr 

@ 

2 

@ 

17 

0 to  @ 

29 

© 

50 

©,  b to  @ 

Precipitation 


some  RW- 


S- 

some  RW- 


Snow  Depth 

2 om 
<1 
<1 
<1 
<1 
0 
0 
0 
0 
0 
0 

2 cm 
2-5 
0 
0 
0 


late  R- 


TABLE  4.  WEATHER  - MARCH  1975 


Date  S 


(Type  I Turbulence) 

4 Mar  9,  brief  ©,@,0 

7 0,  @ 

8 9,0 

9 9,  0 to  ©,  @ 

10  @ 


©,  0 to  @ 


# to  0 


0 to  9 , 0 


(Type  II  Turbulence) 


5 Mar 

6 


9,  0,  @ 


©to  @ 


(Type  III  Turbulence) 


14  Mar 

15 


© to  0 to  @ 


@ to  0 
0 to  i 
@ to  0 


9 to  @ 


9 to 


some  ©,9 


Precipitation 


some  SW- 


SW-,  R- 


some  SW- 


SW-,  s,  s- 


late  RW- 


L,  P,  SW- 


SW-,  S- 


ZRW- 


R- , L- 
late  R- 
P,  brief  R W- 
late  ZR-,  P 


Snow  Depth 


2 cm 


<1-7  cm 
4 
<1 


<1-7  cm 
7-5 


TABLE  6.  WEATHER  - MAY  1975 


Date 

Sky 

Precipitation 

(Type  I Turbulence) 

2 May 

0 to  0,9,0,® 

3 

9 , © to  © 

4 

® 

early  R- 

5 

© to  0 to  9,® 

9 

9 to  0 to  9 

10 

9 to  0 to  9 

11 

9to  @ 

14 

0,®  to  9 to  @ 

16 

0,  @ to  9 

H,  E 

19 

9 to  0 to  @ 

H 

21 

© 

(Type  II  Turbulence) 

7 May 

9 to  0 

13 

0 to  0 to  @ 

some  T 

22 

9,  0,  ®,  0 

F,  H 

23 

© to  0 to  0 

F,  K 

24 

0,some  0,9,® 

early  TRW- 

25 

9,  0 

27 

@ to  OjSorae  0 

28 

0 to  9 

29 

0,  9 

late  RW- 

(Type  III  Turbulence) 

1 May 

® 

RW-,  F 

6 

@ 

RW-,  R-,  F 

12 

@ 

RW-,  F 

17 

© to  ® 

18 

@ 

F,  H 

(continued) 


Snow  Depth 
0 all  month 
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TABLE  6.  (continued) 


Bate 

Siy 

Precipitation 

Snow  Bepth 

20  May 

@ to  ©to  @ 

H,  E 

0 all  month. 

21 

0 

H 

26 

0 

late  RW- 

30 

0 

F,  some  TRW- 

31 

@ 

TRW-,  RW-,  H,  F 
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TABLE  7 


STATISTICS  FOR  LOG  C^2  (TYPE  I TURBULENCE) 


Month. 

Time 

# Pts 

Mean 

Feb 

night 

72 

-2.1293 

day 

59 

-0.7599 

dusk 

28 

-1.5885 

dawn 

32 

-2.2593 

Mar 

night 

157 

-2.3290 

day 

109 

-0.7031 

dusk 

63 

-1.8647 

dawn 

56 

-2.3422 

Apr 

night 

180 

-1.6654 

day 

150 

-0.7913 

dusk 

70 

-1.5856 

dawn 

77 

-1 .7922 

May 

night 

198 

-1.4911 

day 

138 

-0.5602 

dusk 

77 

-1.5285 

dawn 

85 

-1.6567 

Std  Dev 

Skew 

Kurtosis 

0.5818 

-0.0193 

1 .8721 

0.4071 

-1.4601 

4.7501 

0.4554 

-0.7393 

2.3638 

0.6082 

-0.0373 

2.2389 

0.6027 

-0.2786 

2.1450 

0.5226 

-0.8821 

3.0298 

0.6043 

0.5247 

2.4447 

0.6345 

0.3018 

2.4577 

0.4278 

-0.3780 

2.9456 

0.6151 

-0.8924 

3.1236 

0.5569 

0.1467 

2.0468 

0.5187 

0.2501 

3.1293 

0.5369 

-0.1877 

2.3367 

0.3703 

-1.5192 

5.9714 

0.4(534 

-0.2176 

2.9851 

0.5986 

0.0812 

2.5258 

TABLE  9.  STATISTICS  POR  LOG  C.Z  ("ALL  WEATHER") 
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Month 

Time 

# Pts 

Mean 

Std  Dev 

Skew 

Kurtosis 

Peh 

night 

303 

-2.1218 

0.6750 

0.0603 

2.3967 

day 

244 

-1.6349 

0.7631 

0.2102 

2.1841 

dusk 

117 

-1.8117 

0.5895 

-0.2495 

2.6038 

dawn 

136 

-2.2997 

0.7113 

-0.0420 

2.3398 

Mar 

night 

425 

-2.0482 

0.7496 

0.1822 

3.1892 

day 

294 

-1.1210 

0.8219 

-0.3946 

1.9218 

dusk 

160 

-1.7303 

0.6911 

0.0951 

2.2917 

dawn 

163 

-2.1935 

0.7008 

0.2681 

3.3321 

Apr 

night 

249 

-1.6394 

0.5079 

-0.2778 

3.3083 

day 

206 

-0.8155 

0.6286 

-0.8228 

2.9301 

dusk 

98 

-1.6653 

0.6370 

-0.0915 

2.4126 

dawn 

101 

-1.7850 

0.5109 

0.1018 

3.0319 

May 

night 

501 

-1.5387 

0.5728 

-0.6395 

3.1334 

day 

373 

-0.8524 

0.6315 

-1.1269 

3.6039 

dusk 

197 

-1.5110 

0.5412 

-0.3310 

3.0026 

dawn 

199 

-1.7145 

0.6337 

-0.2969 

3.0355 

---  -- 


duslc 


Type  II 


"All  Y/eather" 


dawn 

all 

night 

day 

dusk 

dawn 

all 

night 

day 

dusk 


0,2227 

0,0404 

0.4576 

0.1630 

0.4582 

0.0338 

0.1362 

0.1977 

0.1307 

0.5367 

0.0473 

0.0430 


0.4055 

0.2787 

1.0682 

0.3202 

0.8739 

0.7141 

0.2361 

0.5361 

0.1282 

0.5494 

0.3421 


dawn 


0.1371 
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APPENDIX  A - AVERAGING  TIMES 

I.  THE  EFFECT  OF  AVERAGING  TIMES  ON  THE  TURBULENCE  PROFILE 

Late  in  the  experiment  a question  arose  concerning  averaging  times. 

The  question  originally  came  up  with  respect  to  unexpected  results  in 
the  probability  distributions,  but  it  was  also  unknown  what  effect  the 
long  averaging  time  had  on  the  characteristics  of  the  diurnal  profile, 
especially  on  the  magnitude  of  the  sunset  dip.  The  30-minute  time  period 
for  evaluating  turbulence  levels  from  raw  data  had  been  specifically 
chosen  to  "smooth  out"  the  effects  of  turbulence  intermittency. 

The  48-hour  period  beginning  1400  EST,  9 May  1975,  was  studied  to 
determine  the  effect  of  averaging  time  on  profile  shape.  This  particular 
period  was  picked  because  it  was  the  longest  segment  of  raw  data  collected 
under  good  weather  conditions  still  available  at  the  time  of  this  analysis. 
Good  weather  conditions  were  required  because  they  gave  the  "best"  tur- 
bulence profiles.  The  strong  peaks  in  the  nighttime  turbulence  were  an 
interesting  byproduct  of  the  choice,  but  not  a determining  factor.  The 
raw  data  were  reduced  to  5-minute  periods  in  addition  to  the  usual  30- 
minute  periods  (see  Figures  Ai  and  A2). 

Comparing  Figures  Al  and  A2,  one  immediately  sees  that  the  longer 

averaging  period  effectively  low-pass  filters  the  data.  The  result  is 

helpful  for  observing  trends,  but  it  appears  that  the  sharper  spikes  and 

dips  in  turbulence  are  lost.  The  longer  time  averages  tend  to  follow  the 

spikes  more  closely  than  the  dips  of  the  shorter  time  periods.  The  item 

of  interest,  the  sunset  dip,  changes  substantially  between  the  two  averaging 

periods^  For  5-minute  periods,  its  values  are  1.6x10"^  and  2.0x10'^ 

“C^m  ^ (for  9 May  and  10  May,  respectively),  versus  30-minute  values 
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of  4.8xlO"J  and  1.25x10  °C  m . The  daytime  turbulence  values  are 


.2  -2/3 
or  m ' 


not  affected  as  much,  but  several  spikes  stronger  than  1.0°C  m are 
"lost"  in  the  longer  averages. 

II.  THE  EFFECT  OF  AVERAGING  TIMES  ON  PROBABILITY  DISTRIBUTIONS 

The  probability  distributions  calculated  from  30-minute  averages  of 
raw  data  did  not  exhibit  the  lognormal  shape  predicted  for  daytime  Type  I 
turbulence.  Three  possible  reasons  were  hypothesized,  of  which  one  was  the 
length  of  the  averaging  period.  To  examine  this  possibility,  a 70-hour 
period  of  data,  collected  under  good  weather  conditions  beginning  at  1400 
EST  on  9 May  1975,  was  analyzed. 

The  raw  data  were  reduced  to  5-periods,  from  which  the  usual  probabil- 
ity distributions  were  calculated.  The  statistics  for  log  c|  are  given  below: 

mean  std  skewness  kurtosis  # ots 


.5331 
1.5803 
.6161 
.5624 

The  numbers  suggest  that  the  daytime  distribution  was  significantly  improved 
by  using  5-minute  averages,  whereas  the  nighttime  distribution  flattened, 
almost  to  the  point  of  becoming  a perfect  uniform  distribution  (compare 
Figures  A3  and  A4  with  CIO  and  Cll).  The  cumulative  probability  distri- 
butions of  the  840  consecutive  5-minute  values  and  the  month  of  30-minute 
values  have  been  plotted  on  a probability  scale  in  Figures  A5-A8.  Figure 
A5  has  the  characteristic  shape  of  a uniform  distribution.  A comparison  of 
Figures  A5  and  A8  shows  that  the  distribution  for  the  5-minute  values  lies 
much  closer  to  a straight  line  (a  Gaussian  distribution).  It  appears,  then, 


that  the  long  30-minute  periods  do  limit  the  ability  to  "observe"  strong 
turbulence. 
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Figure  A4.  Discrete  probability  (+)  and  cumulative  probability  (solid  1 
turbulence  for  9-11  May  1975  (day  values):  5-minute  averages 


i>> 

p 

•H 

H 

•H 

& 

o 

(4 

ft 


01 

<i) 

3 


P 

■a 

•H 

P) 


in 

o- 

crv 


i 

cr» 


u 

o 

<H 


Pi 

O 

•H 

P 

O 

rt 

3 


t>> 

P 

•H 

H 

•H 

.O 


m 


d o 


& 

o 


SP 


P U 
ft  <u 
> 


<D  ci 
i> 


•p  d) 


p p 
uS  3 


iH 

3 


•H 

a 

3 I 

o in 


in 


a) 

P 


a 


•H 

Pm 


86 


Figure  A6.  Cumulative  probability  function  for  9-11  May  1975  (day  values): 
5-minute  averages. 
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Figure  A8.  Cumulative  probability  function  for  May  1975,  Type  I turbulence  (day 
values):  30-minute  averages. 


APPENDIX  B - SENSOR  NORMALIZATION  ERRORS 


; i 


\ 


t 


i 

k 
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The  length  of  the  thin  platinum  wire  which  is  used  as  the  temperature 
sensing  element  of  microtemperature  probes  varies  slightly  from  unit  to 
unit.  This  causes  the  probes  to  have  different  gain  factors,  which  results 
in  nonuniform  voltage-to-temperature  conversions.  The  resistance  of  the 
platinum  wire  is  measured  and  compared  to  a standard  value  for  normalization. 
If  the  resistance  measurements  are  in  error  for  any  reason,  the  probes  will 
not  be  correctly  normalized.  It  was  feared  that  this  could  grossly  affect 
the  calculation  of  C2.  The  following  analysis  explores  the  problem. 

Since  the  probe  voltage  is  AC  coupled  to  ambient  temperature  changes 
(the  timr  constant  is  approximately  100  sec), 

< voltage  >:0  (Bl) 

or,  since  the  temperature  T is  directly  proportional  to  voltage, 

< T > = 0 (B2) 

Then 

aj  = < (Tj)2>  , (B3) 

(a2')2  = < (T  1 )2  > , (B4) 

2 

and 

C ' = < (T  ) (T  ' ) > (B5) 

12  12 

where  the  primes  indicate  values  affected  by  normalization  errors.  We 
can  write 

(T  1 )2  = (T  )2  (1+c)  (B6) 

2 2 

where  e accounts  for  possible  differences  between  probe  calibrations.  C£ 


t 

I 


\ 

i 

1 

, 

<’l 

! 
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is  related  to  the  temperature  variances  and  covariances  by 

Ci  « 02  + (0  1 ) - 2C  1 (B7) 

t 1 2 12 

From  B4-B6  we  find  that 

(a;)2  = 2 (1+s)  (B8) 

and 

C,a'  * C12(1+E)I/2  (B9) 

Substituting  back  into  B7,  we  have 

C?  « a2  + a 2 (1+e)  - 2C  (1+e)1^2  (BIO) 

t 1 2 12 

If  it  is  assumed  that  e is  "sinal  1",  then 

(1+e) 1/2  - l+(e/2)  (BIT) 

Note  that  for  e as  large  as  0.5,  the  approximation  is  less  than  3%  in  error. 
Using  (Bll),  (BIO)  becomes 

C2  « o 2 + 022  (1+e)  - 2C  (1+4-  e) 

=a  2 + a 2 - 2C  + e(0  2-  C ) (B12) 

1212212 

The  error  has  now  been  isolated  to  a single  term.  Suppose  the  desired 
accuracy  is  only  25%;  then 

e(o  2 - C ) < (0  2 + 0 2 - 2C  ) /4 
2 12  1 2 12 

or  2 

0 + 0 2 - 2C 

e < 1 -i L2 

4(0  2 - C ) (B13) 

2 1 2 

From  turbulence  theory,  we  expect  o^2  » o*  for  sufficient  averaging  times. 
Equation  (B13)  then  becomes 


2 


= 2(a  2 - C ) 
1 12 


4(o  2 - C ) 
1 12 


= 0.5  (B14) 

Under  the  assumptions  which  have  been  made,  errors  in  sensor  cali- 
bration up  to  50  percent  will  cause  less  than  a 25  percent  error  in  the 
calculation  of  Ct-.  Such  a gross  error  in  normalization  is  unlikely;  a 10 
percent  error  is  more  plausible.  This  would  result  in  an  error  in  Ct2  of 
only  5 percent.  Therefore,  normalization  errors  should  present  no  problem 
to  data  analysis. 


APPENDIX  C - PROBABILITY  DISTRIBUTIONS 
This  appendix  contains  the  plots  of  the  probability  distributions. 
The  distributions  were  calculated  for  each  turbulence  type  from  the  30- 
minute  Cf  values,  divided  into  night,  day,  dusk  and  dawn.  The  cumulative 
probability  functior  (CPF)  is  plotted  as  a continuous  line,  while  the 
discrete  probability  function  (OPF)  is  plotted  as  a sequence  of  points 
which  give  the  probability  of  a value  occurring  within  the  discrete 
interval  around  that  point.  Another  common  function,  the  probability 
density  function,  can  be  easily  calculated  from  the  discrete  probability 
function  by  dividing  the  probabilities  by  the  interval  (width)  that  each 
point  represents. 
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Probability 


Figure  C2.  Discrete  probability  (+)  and  cumulative  probability  (solid  line)  for 
turbulence:  February  1975,  day. 


Figure  C6.  Cumulative  probability  functions  for  Type  I turbulenc 
dusk  and  dawn. 


Figure  G8.  Discrete  probability  (+)  and  cumulative  probability  (solid  line)  for  Type 
turbulence:  April  1975,  day. 


T 


Figure  Oil.  Discrete  probability  (+)  and  cumulative  probability  (solid  line)  for  Type 
turbulence:  May  1975,  day. 


Figure  Cl  3.  Discrete  probability  (+)  and  cumulative  probability  (solid 
Type  II  turbulence:  February  1975,  nigbt. 


Figure  Cl 4.  Discrete  probability  (+)  and  cumulative  probability  (solid  line)  for 
Type  II  turbulence:  February  1975,  day. 
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Figure  Cl  7.  Discrete  probability  (+)  and  cumulative  probability  (solid  line)  for 
Type  II  turbulence:  March  1975,  day. 


.Figure  Cl 8.  Cumulative  probability  functions  for  Type  II  turbulence:  March  1975 
dusk  and  dawn. 


probability  (+)  -and  cumulatire  pr obab 
turbulence:  April  1975,  night. 


Figure  C20.  Discrete  probability  (+)  and  cumulative  probability  (solid  line)  for 
Type  II  turbulence:  April  1975,  day. 
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Figure  C25.  Discrete  probability  (+)  and  cumulative  probability  (solid 
turbulence  in  "all  weather":  February  1975,  night. 
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Figure  C29.  Discrete  probability  (+)  and  cumulative  probability  (solid  line)  for 
turbulence  in  "all  weather":  March  1975,  day. 


Figure  C30.  Cumulative  probability  functions  for  turbulence  in  nail  weather” 
March  1 975 , dusk  and  dawn. 


Figure  C34.  Discrete  probability  (-f)  and  cumulative  probability  (solid  line)  for 
turbulence  in  "all  weather":  May  1975,  night. 


dawn 


Figure  031.  Discrete  probability  (+)  and  cumulative  probability 
turbulence  in  "all  vre  airier" : April  1375,  night. 


Figure  C32.  Discrete  probability  (+)  and  cumulative  probability  (solid 
turbulence  in  "all  weatber” : April  1975,  day. 


Figure  C35.  Discrete  probability  (+)  and  cumulative  probability  (solid  line)  for 
turbulence  in  "all  weather”:  May  1975,  day. 


MISSION 
of 

Rome  Air  Development  Center 


RADC  plans  and  conducts  research,  exploratory  and  advanced 
development  program  in  command,  control,  and  communications 
(C3)  activities,  and  in  the  C3  areas  of  information  sciences 
and  intelligence.  The  principal  technical  mission  areas 
are  communications,  electromagnetic  guidance  and  control, 
surveillance  of  ground  and  aerospace  objects,  intelligence 
data  collection  and  handling,  information  system  technology, 
ionospheric  propagation,  solid  state  sciences,  microwave 
physics  inid  electronic  reliability,  maintainability  and 
compatibility. 


